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Abstract 


We study the nonequilibrium dynamics leading to the formation of topologi¬ 
cal defects in a symmetry-breaking phase transition of a quantum scalar field 
with self-interaction in a spatially flat, radiation-dominated Friedmann- 
Robertson Walker Universe. The quantum field is initially in a finite- 
temperature symmetry-restored state and the phase transition develops as 
the Universe expands and cools. We present a first-principles, microscopic 
approach in which the nonperturbative, nonequilibrium dynamics of the quan¬ 
tum field is derived from the two-loop, two-particle-irreducible closed-time- 
path effective action. We numerically solve the dynamical equations for the 
two-point function and we identify signatures of topological defects in the 
infrared portion of the momentum-space power spectrum. We find that the 
density of topological defects formed after the phase transition scales as a 
power law with the expansion rate of the universe. We calculate the equi¬ 
librium critical exponents of the correlation length and relaxation time for 
this model and show that the power law exponent of the defect density, for 
both overdamped and underdamped evolution, is in good agreement with the 
“freeze-out” scenario proposed by Zurek. We also introduce an analytic dy¬ 
namical model, valid near the critical point, that exhibits the same power-law 
scaling of the defect density with the quench rate. By incorporating the realis¬ 
tic quench of the expanding Universe our approach illuminates the dynamical 
mechanisms important for topological defect formation and provides a pre¬ 
liminary step towards a complete and rigorous picture of defect formation 
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in a second-order phase transition of a quantum field. The observed power 
law scaling of the defect density with the quench rate, calculated here in a 
quantum field theory context, provides evidence for the “freeze-out” scenario 
in three spatial dimensions. 
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I. INTRODUCTION 


A. Motivation 

Phase transitions in the early Universe played a decisive role in shaping the Universe 
we observe today. In the hot early Universe it is likely that the broken symmetries of 
the standard model of particle physics were partially or totally restored in a grand unified 
theory (GUT) and that a number of phase transitions occurred as the Universe expanded 
and cooled. These include transitions at the GUT scale T^. k, 10^“^ — 10^® GeV, at the 
electroweak scale 10^ GeV and the color deconhnement and/or chiral phase transitions 
of QGD with Tc 100 MeV. An inflationary phase transition, a possible explanation for 
the flatness, horizon and monopole problems of the standard Friedmann-Robertson-Walker 
(FRW) cosmology, may have occurred at the GUT scale or before. At still earlier epochs 
near the Planck energy, ~ 10^® GeV, candidate theories of quantum gravity may allow for 
a phase transition to occur and give rise to the classical properties of spacetime described 
by general relativity. These first moments of the Universe contain clues to the solution 
of such outstanding problems as the origin of large-scale structure and matter-antimatter 
asymmetry and the nature of dark matter. The study of the consequences of cosmological 
phase transitions may therefore hold the key to our understanding of the Universe. 

In a held theory with degenerate vacua, a symmetry-breaking phase transition can pro¬ 
duce topological defects. Relics of the high-temperature symmetric phase of the theory, 
topological defects are topologically stable held conhgurations that are locally trapped in an 
excited state above the vacuum. Examples of defects include vortices in type II supercon¬ 
ductors, superhuid IIehum-3 and IIeIium-4, and cosmic strings formed in the early Universe. 
Topological defects are classihed by the homotopy groups n„(M) of the vacuum manifold, 
M. In three spatial dimensions simple defects are domain walls if no(M) ^ 1, strings if 
ni(M) 7 ^ 1, and monopoles if 112 (M) 7^ 1 |]1|1- 

Kibble was the hrst to show that topological defects are a generic feature of cosmological 
phase transitions [0. In the course of a symmetry-breaking phase transition, as the held 
decays to the stable vacuum, it can choose the same vacuum state on length scales only as 
large as the correlation length. In the laboratory, if the phase transition proceeds slowly, 
the correlation length is bounded only by the size of the system and the held will ehectively 
choose a homogeneous vacuum state. In the early Universe, however, the correlation length 
is bounded by the size of the particle horizon. With correlated regions limited in size by 
causality. Kibble argued that phase transitions in the early Universe necessarily leave a 
domain structure of vacuum states and form topological defects. 

Topological defects formed in the early Universe can have profound consequences on the 
subsequent spacetime evolution. In the standard FRW cosmology, GUT scale monopoles, 
even if formed with an initial density equal to Kibble’s lower bound of one defect per horizon 
volume, contribute more than 10^^ times the largest possible density of the Universe consis¬ 
tent with observations. The overproduction of monopoles is a puzzle within the standard 
FRW cosmology and provided part of the original motivation for an early epoch of inflation 
0. Other topological defects such as cosmic strings are viable candidates for the seeds 
of structure formation and would produce a characteristic signature in the angular power 
spectrum of the cosmic microwave background. Satellites such as MAP and Plank designed 
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specifically to measure fluctuations in the cosmic microwave background radiation to high 
accuracy are scheduled to be launched in the near future. The data obtained from these 
missions is expected to clarify the precise role of topological defects in the origin of structure 
formation in the Universe. 

The importance of topological defects as candidates for the seeds of structure formation 
provides strong motivation to extract predictions from defect models that can be compared 
with observations. These predictions generally contain three main ingredients: the den¬ 
sity of topological defects immediately following the completion of the phase transition, the 
subsequent evolution of the defect distribution and the calculation of physical observables 
such as cosmic microwave background radiation anisotropies and polarization. While both 
numerical and analytical work has been done on the evolution of the defect distribution and 
the extraction of physical predictions, less attention has been focused on quantitative pre¬ 
dictions of the initial defect density which provides the starting conditions for the evolution 
of the defect distribution. 

The dynamic origins of defect formation are also themselves of considerable physical 
interest. Formed in out-of-equilibrium processes, defects offer insight into the nature of crit¬ 
ical dynamics. In addition, the ubiquitous nature of phase transitions links the study of the 
early Universe with the laboratory and allows the observation of cosmologically important 
mechanisms in condensed matter systems Q |^. Recently, researchers have exploited these 
similarities by observing defect formation in liquid crystals 0 and superfluid Helium [0. 
In either condensed matter systems or the early Universe, the density of topological defects 
immediately following the phase transition is of substantial physical importance. 


B. The Kibble and Zurek Mechanisms of Defect Formation 

The hrst estimate of the initial defect density in a cosmological context was made by 
Kibble p|. The basic ingredients of the Kibble mechanism are causality and the Ginzburg 
temperature, Tq- The Ginzburg temperature is dehned as the temperature at which thermal 
fluctuations contain just enough energy for correlated regions of the held to overcome the 
potential energy barrier between inequivalent vacua, 

hTG^iiTcf^FiTc), ( 1 . 1 ) 

where AF is the difference in free energy density between the true and false vacua and ^ is 
the equilibrium correlation length. In the Kibble mechanism, the length scale characterizing 
the initial defect network is set by the equilibrium correlation length of the held, evalu¬ 
ated at Tg- In a recent series of experiments 0, the Kibble mechanism was tested in the 
laboratory. The results, while conhrming the production of defects in a symmetry-breaking 
phase transition, indicate that ^{Tg) does not set the characteristic length scale of the initial 
defect distribution. These experiments were suggested by Zurek, who criticized Kibble’s use 
of equilibrium arguments and the Ginzburg temperature, Tg 0. 

Gombining equilibrium and non equilibrium ingredients, Zurek ohers a “freeze-out” pro¬ 
posal to estimate the initial density of defects. In the Zurek proposal, above the critical 
temperature, the held starts oh in thermal equilibrium with a heat bath. As the temper¬ 
ature of the bath is lowered adiabatically, the held remains in local thermal equilibrium 
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with the heat bath. Near the phase transition, the equilibrium correlation length and the 
equilibrium relaxation time of the held grow without bound as 


e = eo|er^ (1.2) 

r = ro|e|“'', (1.3) 


where e characterizes the proximity to the critical temperature, 

( 1 - 

and /i and v are critical exponents appropriate for the theory under consideration. The 
quench is assumed to occur linearly in time 



so that for f < 0, the temperature of the heat bath is above the critical temperature and the 
critical temperature is reached at t = 0. The divergence of the equilibrium relaxation time as 
the heat bath approaches the critical temperature is known as critical slowing down. Critical 
slowing down results from the hnite speed of propagation of perturbations of the order 
parameter. As the correlation length diverges, small perturbations of the order parameter 
lowering of the temperature) take longer to propagate over correlated regions and 
therefore it takes longer to reach equilibrium. As the critical temperature is approached 
from above there comes a time | t* \ during the quench when the time remaining before the 
transition equals the equilibrium relaxation time 


t* 1 = T{t*). 


( 1 . 6 ) 


Beyond this point the correlation length can no longer adjust fast enough to follow the chang¬ 
ing temperature of the bath. At time t* the dynamics of the correlation length “freezes”. 
The correlation length remains frozen until a time t* after the critical temperature is reached. 
In Zurek’s proposal the correlation length at the “freeze-out” time f^sets the characteristic 
length scale for the initial defect network. Solving for the value of the correlation length 
at t*, the frozen correlation length, and therefore the initial defect density, scale with the 
quench rate as 

^(t*)r^TQ^. (1.7) 

Aspects of this scenario have been verihed using 1-dimensional and 2-dimensional simula¬ 
tions of phenomenological time-dependent Landau-Ginzburg equations |^. Preliminary 
experimental efforts to test the Zurek prediction in condensed matter systems are incon¬ 
clusive, lacking reliable error estimates and a mechanism to vary the quench time scale 

i- 

The Kibble prediction for the initial length scale of the defect distribution is inconsistent 
with experiment, yet the alternative “freeze-out” proposal, raises many conceptual questions. 
The “freeze-out” scenario of Zurek relies on phenomenological ideas derived from experience 
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with classical Landau-Ginzburg systems. It is not clear how much of this picture, if any, is 
applicable to quantum helds which undergo phase transitions in the early Universe. Quantum 
helds bring new conceptual issues. In particular, consideration of decoherence is necessary 
to understand when and how the quantum system at the late stages of the transition can be 
described by an ensemble of classical defect conhgurations. Even in the classical context, the 
use of equilibrium critical scaling in a fully dynamical setting is an approximation. Critical 
dynamics is generally much richer and less universal than equilibrium critical behavior p3 


For a system coupled to a rapidly changing environment, the correlation length evolves 
through the dynamical equations of motion of the held. Worse still, the relaxation time is 
not always well-dehned. In addition, a realistic bath consists of interactions between the 
system and other helds in the Universe. Physical interactions between the bath and the 
system produce interesting and nonequilibrium behavior such as noise and dissipation []n 


The reduction of the complicated bath-system interaction to one parameter, the quench time 
scale, and a prescribed linear time-dependence in the ehective mass of the time-dependent 
Landau Ginzburg equation neglects these potentially important processes. 


C. Foundational Issues 

Despite signihcant ehort, the physical mechanisms important for the formation of topo¬ 
logical defects in the early Universe and condensed matter systems are not well understood. 
It is useful to explore three interwoven and complementary themes that are relevant to this 
problem. These are (i) equilibrium vs. nonequilibrium, (ii) microscopic vs. macroscopic, 
and (iii) quantum vs. classical. 


1. Equilibrium vs. Dynamical 

The Zurek and Kibble mechanisms both rely on aspects of equilibrium physics. However, 
as emphasized by Zurek, the formation of topological defects is a nonequilibrium process. If 
aspects of the “freeze-out” scenario are correct then it must be viewed as an approximation to 
a more fundamental and as yet unexplored understanding. The initial density of topological 
defects formed in far-from-equilibrium phase transitions may be determined by very different 
mechanisms. 


2. Microscopic vs. Macroscopic 

Topological defects appear as conhgurations of the classical order parameter held in 
phenomenological Landau-Ginzburg theories and in the classical limit of symmetry-broken 
quantum held theories. A hrst-principles approach to the formation of topological defects in a 
nonequilibrium quantum held phase transition requires an understanding of the microscopic 
origin of macroscopic dynamical critical behavior like critical slowing down. Macroscopic 
domains of correlated vacuum must be identihed from the microscopic quantum held. Ther- 
malization occurs through microphysical couplings between the quantum held composing 
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the system and other helds composing an environment and mnst be nnderstood to prodnce 
a physical bath-system interaction. 


3. Quantum vs. Classical 


Phase transitions occnr in qnantnm held theory and in condensed matter systems. In a 
second-order phase transition, the correlation length diverges at the critical temperature and 
the dynamics is dominated by low frequency modes. In equilibrium, modes with uj « T 
have essentially a Maxwell-Boltzmann distribution and it might be expected that second- 
order phase transitions are dominated by classical behavior. In a dynamical setting an exact 
equilibrium distribution does not exist and a necessary indication of classical behavior is the 
emergence of a positive dehnite probability distribution function from the density matrix 
of the quantum held 0- In quantum critical systems decoherence may modify critical 
scaling exponents and change the coarsening of domains from that of classical systems. 
Even the concepts used to describe phase transitions in condensed matter may not apply in 
the quantum context. For example, the dehnition of a topological defect in a quantum held 
theory with arbitrary quantum state is unknown. 


Phase transitions and topological defects occur in both classical and quantum systems, 
under equilibrium and nonequilibrium conditions and are depicted by phenomenological 
and microscopic theories. They therefore provide a theoretical vantage point from which 
to further elucidate the connections between these complementary and disparate concepts. 
The formation of topological defects is an arena in which these general yet basic issues of 
physics may be studied in a concrete and useful manner. 


D. Quantum Field Dynamics and the Formation of Topological Defects in the Early 

Universe 


A complete understanding of the physical issues involved in the formation of topological 
defects in a second-order phase transition in the early Universe requires a hrst-principle 
approach to the nonequilibrium dynamics of quantum fields, a realistic treatment of the 
interaction of the quantum held both with gravity and with other helds that constitute an 
environment, and the identihcation of classical defect conhgurations from the quantum held 
system. 

This work will concentrate on the dynamics of the phase transition important for defect 
formation. While symmetry restoration in hnite-temperature quantum held theory has been 
known since the early work of Kirzhnitz and Linde |0 , most previous ehorts have focused 
on the equilibrium aspects of the transition. Techniques such as the hnite-temperature 


ehective potential JTl 0 and the renormalization group |T^ have been developed to deduce 
equilibrium critical properties such as the order of the phase transition and its critical 
temperature. However, equilibrium techniques are inadequate to study the dynamics of 
the phase transition. The use of equations of motion for the mean held derived from the 
hnite-temperature ehective potential was criticized in |]^. In general the use of equations 


generated from the hnite-temperature ehective potential in a dynamical setting results in 
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unphysical solutions. Although equilibrium techniques are clearly inappropriate, solving 
the full equations of motion for an interacting field theory is generally impossible, even 
numerically. This difficulty is partially overcome by the development of approximation 
schemes which allow for the evolution of a restricted set of correlation functions of the 
quantum field theory . These methods have been applied to a number of dynamical 

problems in quantum field theory including post-inflationary reheating |]T^ (and references 
therein) and spinodal decomposition []^ |T^ . 

The problem of defect formation in a nonequilibrium second-order phase transition of a 
quantum field has also recently received attention [^] . The results are promising but 

the studies are incomplete. In these previous approaches the phase transition is incorporated 
through an ad hoc time dependence of the effective mass of a free field theory: an instability 
in the theory is induced when the mass becomes tachyonic. The use of a prescribed time 
dependence of the effective mass, while providing a convenient analytic model, lacks physical 
justification. The neglect of interactions confines the applicability of these approaches to 
very early times before the field amplitude grows substantially. They are therefore unable to 
account for the back reaction which is necessary to stabilize domain growth and shut off the 
spinodal instabilities of the phase transition. In addition, defects formed during the linear 
stages of the phase transition are transient and not likely to survive to late times. 


E. Outline 

In this paper we analyze the formation of correlated domains of true vacuum in the 
time evolution of a quantum scalar field during a second-order phase transition initiated by 
cooling of a radiation-dominated FRW Universe. In Section II we discuss the formalism used 
to follow the dynamics of a quantum field through the phase transition. We also present the 
model, a derivation of the equations of motion for the two-point function of the theory, and 
discuss renormalization, initial conditions and numerical parameters used in the numerical 
simulation. We then provide a dynamical description of the phase transition using results 
of the numerical simulations. Section III begins the discussion of domains and presents the 
argument that domains are determined by a peak in the Fourier space structure function 
kfG{k, t). Section IV discusses the power law scaling of the size of domains with the quench 
rate. The equilibrium critical exponents of the correlation length and relaxation time are 
calculated in both the underdamped and overdamped cases and the power law scaling of 
the defect density with the quench rate predicted by the “freeze-out” proposal is shown to 
be in good agreement with the numerical simulations. An analytical model valid for slow 
quenches and near the onset of the instability is introduced and the power law exponent in 
the analytic model is found to be the same as the numerical simulations. Section V provides 
a summary and discussion of these results and presents possible directions for further study. 


II. PHASE TRANSITION DYNAMICS 

It is common to model critical dynamics with a classical, phenomenological, time- 
dependent Landau-Ginzburg equation for an order parameter T, 






(2.1) 


s,<i>(s,«) = -r—+ e, 

where F[^] is a phenomenological free energy density for the order parameter, F is a phe¬ 
nomenological dissipative coefficient and is a stochastic term incorporating thermal fluc¬ 
tuations of the environment . 

Even if the order parameter is of quantum origin, as e.g. in the phase of the wavefunction 
for liquid Helium-4, the Landau-Ginzburg equation is rarely derived logically from the un¬ 
derlying quantum dynamics of the system. In condensed matter systems, to compensate for 
insufficient microscopic information, great care with physical intuition goes into choosing the 
order parameter and its equation of motion. In experimentally inaccessible environments, 
such as the early Universe, it is not a priori obvious what the order parameter, or its dy¬ 
namics, should be. In situations where phenomenological approaches are inadequate, it is 
necessary to work with the fundamental quantum dynamics of the helds. 

A hrst-principles approach to the quantum dynamics of phase transitions avoids ad hoc 
assumptions about the dynamics of the correlation length and the effect of the quench. The 
system simply evolves under the true microscopic equations of motion. We can therefore 
explore many details of critical dynamics that are inaccessible in phenomenological theories. 

Solving the exact dynamics of interacting quantum helds expressed, for example, through 
the Heisenberg equations of motion for the held operator is a very complicated problem. 
Instead, we seek a truncation of the full degrees of freedom that allows for an accurate 
modeling of the dynamics of the phase transition over the time scales of interest. 

Consider a classical <I>^ scalar held theory with > 0 and self-coupling A. The potential 
I/[<I)] for the held is 


^ ^ 2 4 

The degenerate true minima of this potential are 


( 2 , 2 ) 


^ min. — i 


m 


(2.3) 


which are nonperturbatively large in the coupling constant. To follow the held as it evolves 
from the unstable vacuum 4) = 0 to the stable one 4) = ±y requires an approach which is 
both non-perturbative in the held amplitude and fully dynamical. 

A useful formalism that is both nonperturbative and dynamical is the two-particle irre¬ 
ducible closed-time-path (2PI-CTP) or in-in ehective action [^. The 2PI-CTP ehective 
action generates real and causal equations of motion for the mean held {^{x,t)) and two- 
point correlation functions {^{x,t)^{y,t')) of the quantum held theory. A diagrammatic 
expansion of the 2PI ehective action, to arbitrary loop order is given in [^]. Calculated to 
all orders in a loop expansion, the 2PI-CTP ehective action and the equations of motion 
derived from it contain all information of the original quantum held theory. However, any 
practical computation requires inclusion of terms only to some hnite loop order which con¬ 
stitutes an approximation |^. In this paper, we neglect three-loop and higher graphs. 


The truncation of the 2PI-CTP ehective action to two-loop order is equivalent to the time 
dependent Hartree-Fock approximation . 
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The equations of motion for the mean held and the two-point function derived from the 
2PI-CTP effective action respect the $ —symmetry of the classical action. Since the 
held starts in a symmetry-restored state above the critical point where 


($) 


initial 


= 0 , 


(2.4) 


the mean held remains identically zero throughout the phase transition. The dynamics of 
the phase transition unfolds through the dynamics of the two-point function. 

The choice to limit attention to the two-point correlation functions in the Hartree-Fock 
approximation limits the time scale over which the evolution is physically reliable. Higher 
order correlations are suppressed in weakly coupled held theories in equilibrium. However, 
second-order phase transitions are characterized by spinodal instabilities that cause corre¬ 
lations to grow. When the correlations grow large enough that the held is sampling the 
stable vacua, dissipative processes are expected to become important and the Hartree-Fock 
approximation will break down 


20 . 


A. The Model 


We consider a scalar held in an FRW spacetime 
breaking classical action 


30]. The held has the symmetry- 


s = J ( 2 . 5 ) 

where g is the determinant of the metric of the classical background spacetime. We assume 
that the stress-energy tensor is dominated by other radiation helds present in the early Uni¬ 
verse. These helds maintain the overall homogeneity and isotropy of the Universe. Small 
deviations in homogeneity and isotropy that may eventually be responsible for the huctu- 
ations observed in the cosmic microwave background radiation are produced in our model 
by the topological defects of the system and appear only at the end of the phase transition. 
We therefore work in the semiclassical test-held approximation (ignoring back reaction of 
the <F held on the spacetime), and we assume that the scale factor has the time-dependence 
of a homogeneous and isotropic, spatially hat, radiation-dominated Universe, 


a(t) 


t + T 


T 


(2.6) 


The expansion of the Universe and the resulting redshifting of the modes act here as a 
physical quench allowing the dynamics of the phase transition to unfold naturally. This is 
in distinction to work which uses an instantaneous change in the sign of the square of the 

g 0. 


mass 


In the Heisenberg representation the held operator can be written as 







-ik-x . 


(2.7) 


where the fk are the complex quantum modes of the held. 
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In a FRW Universe, the two-loop 2PI equation of motion in comoving cosmological time 
for the mode function with comoving momentum k is 


dt‘^ a(t) dt a^(t) 


+ ^G(i, i)) fk{t) = 0. 


(2.8) 


where 

is the equal-time limit of the two-point correlation function and 

fkdi) = cotli(^uii(0)) 


(2.9) 


( 2 . 10 ) 


is a constant factor incorporating thermal initial conditions with temperature T = 1/j3. The 
effective mass of the system is 

t) (2.11) 

and the initial frequency is 

wl{S)) = e + m%^{S)), ( 2 . 12 ) 

where mgjy(O) is the initial hnite-temperature effective mass in the Hartree-Fock approxi¬ 
mation. Since we work in a radiation-dominated FRW universe, the scalar curvature, R, is 
zero and the conformal coupling constant ^ may be ignored. 


1. Renormalization 


The equal-time limit of the two-point function is divergent and must be regularized. 
A simple regularization method, amenable to a numerical simulation, is to implement an 
ultraviolet cutoff in physical spatial momentum. A suitably regularized expression for the 
two-point function must be independent of this cutoff. Our renormalization scheme follows 


that of 1^ 


It is a general feature of quantum held theory that the bare quantities appearing in the 
classical Lagrangian are not observable but are “dressed” by interactions. In the Hartree 
approximation, the effect of interactions is encoded in the self-consistent effective mass. 
Therefore, there are no counterterms to absorb divergences and the equation for the mode 
functions, Eq. (p^.8|) , must be hnite [^. We therefore hx the renormalization scheme with 
the condition 


- = -“/i + 


(2.13) 


The cutoff dependence of the bare variance is obtained by considering a WKB-type solution 
to the mode function equation. Identifying the second order adiabatic mode functions from 
the WKB solution of the mode function equation, the subtracted two-point function is 
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(2.14) 


Gs{t, ^ (hfkC’kiP) - 


1 9{k — k) 2 ^ 

+ 4P 


In combination with the renormalization condition, this subtraction can be implemented by 
a shift in the bare parameters of the theory 


2 , Aj A" 

+ „ „, , = m. 


Ibvr^ a^{t) 


1 + 


Aft 


IbTT^ 


ln(A/ k) 


Aft — 


1 - Te> ln(A/K) 


(2.16) 

(2.16) 


The shift is time independent as long as the cutoff A and the renormalization scale k are 
implemented in terms of the physical momentum 


A = Aoa(t) 

K = Koa{t) 


(2.17) 

(2.18) 


where A and k are comoving and Aq and Kq are physical quantities. The renormalized mode 
function equation is 


df^ 


d{t) d k‘^ 
a{t) dt a'^{t) 




+ t) 


fkif) = 0. 


(2.19) 


In following sections we will drop the renormalization subscripts for clarity and it is to be 
understood that we are working with renormalized quantities. 


2. Initial Conditions 


The quantum held is assumed to be initially in a state of thermal equilibrium. In 
this model with a tachyonic tree-level mass, symmetry is restored by hnite temperature 
corrections. The initial effective mass ^^e solution of the (renormalized) equation 


^eff = 


-m? + 


A 


k'^dk 


both 

2^k'‘ + ml„ 


2k 


+ e{k 


m 


K 


eff 


4fc3 


( 2 . 20 ) 


In the high temperature and small-A limit this yields 

AT2 


m 


eff 


= —m -I- 


24 


( 2 . 21 ) 


which is a result familiar in hnite-temperature held theory 


In an expanding FRW Universe, exact thermal equilibrium will only persist for confor¬ 
mally invariant helds. If the expansion rate is small relative to internal collisional processes 
of the held then there is an approximate notion of equilibrium |]32 
transforming to conformal time, r], dehned by 


This is evidenced by 
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dt = a{ri)dr], 


( 2 . 22 ) 


and performing a mode redefinition 

fkiv) = fkivMv) 

The conformal mode fnnction eqnation is now 

(^ + fc^ + «^(h)("i^ + ^G'(p,r7))) fk{T]) = 0. 


(2.23) 


(2.24) 


If the Universe is slowly expanding, a WKB-type solntion is appropriate and a low- 
adiabaticity trnncation of the instantaneons WKB freqnency is snfficient. The zeroth- 
adiabatic order solntion to eqnation (|2.24|) is given by 

fiv) = 


^ 3 -* P Wkdri' 


pMk 
= k^ + 

This leads to the following initial conditions for the mode fnnctions in cosmic time t 


wl = k'^+ a{rifmljj. 


(2.25) 

(2.26) 


fkiO) = 


A-(o) = , } ^ , 

a(0)y2M;fc(0) 

(2.27) 

■jo) "“‘WAW- 

(2.28) 


Slow expansion assnmes that the natnral freqnency of the kth mode is faster than the 
expansion rate of the Universe, 


w^fc(O) > (2.29) 

2r 

The adiabatic eqnilibrinm approximation wonld fail for the lowest k modes. However with 
high temperatnre initial conditions the low k modes are not a dominant part of the spectrnm. 


3. Numerical Parameters 

The model described by Eq. (|27^) and ( |2.6| ) is characterized by 6 parameters, 

(m, r. A, T, A, k). (2.30) 

We choose nnits of energy in which 

rP = 1.0. (2.31) 

The initial rate of expansion or the initial Hnbble constant is controlled by r. 
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(2.32) 


H{0) 


a(0) 


1 


In the simulations we conducted the range of the r parameter was 


0.01 < r < 100 


(2.33) 


and the self-coupling is strong, 


A = 0.1. (2.34) 

The initial temperature T is chosen so that the value of the initial effective mass is 

of order unity. More specifically we choose 

m|//(0) = -1-0 + = 607 (2.35) 

T = 20.0 (2.36) 

The values of the initial effective mass the coupling A and expansion rate r, were 

chosen so that the simulations of the phase transition completed on numerically accessible 
timescales. Due to the renormalization scheme, both the comoving cutoff A and the comoving 
renormalization scale k increase with the scale factor. The initial value 


Ao = 340 (2.37) 

was chosen as the lowest value that maintained cutoff independence of the initial effective 
mass. Insensitivity to the cut-off in the dynamical simulations was verihed by doubling Aq 
and observing no change in the output plots of the time-dependent effective mass. The 
initial value of the renormalization scale 


Ko = 1.0 (2.38) 

was chosen so that the renormalization scale was always above the maximum momentum 
of the unstable modes. The coupled, nonlinear system of mode function equations with 
the given initial conditions was solved numerically using an adaptive stepsize, fifth-order 
Runge-Kutta code. Mode integrals were performed using a simple simpson rule with a 
uniform momentum binning 


h — flkhin, 

h- - — 

^bin — Y ’ 
-^0 


(2.39) 

(2.40) 


where Lq = 100.0 is the effective size of the system and n is the total number of modes. 
Insensitivity to the momentum binning was verihed by reducing kun and observing no change 
in output. The number of modes n varied from 10^ to 10®. Run times varied from hours to 
days on a DEC 500 MHz workstation which corresponds to dynamical time scales of t = 2 
to t = 100. 
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4- Analysis 


The results of a typical simulation are shown in Figs. 1 and 2. In Fig. 1 the renormalized 
effective mass is shown as a function of cosmological time for quench parameter r = 1.0. 
The phase transition begins when hrst becomes negative. When is negative, 

modes with physical momentum ^ < Me// have imaginary frequencies and begin to grow. 
This indicates the onset of the spinodal instability which is characteristic of a second-order 
phase transition. In the early stages of the phase transition, the evolution of the effective 
mass is dominated by the redshifting of stable modes and the effective mass decreases. As 
the phase transition proceeds, more modes redshift into the unstable momentum band and 
the amplitude of unstable modes continues to grow. Eventually, the redshift of stable modes 
balances the growth of unstable modes and the effective mass increases. As the effective 
mass passes through zero from below, the field reaches the spinodal point. Beyond the 
spinodal point all modes are stable. As the effective mass continues to grow, non-linear 
thermal and other collisional processes are expected to be important and the Hartree-Fock 
approximation of the dynamics of the two-point correlation function breaks down (see for 


example ||I^). 

In Fig. 2 we show the Fourier space structure factor 


S{k,t) = k^Gkit.t) 


(2.41) 


at various times during the phase transition. As argued in the next section, defects and 
domains can be identified in the low-k structure of S{k,t). As the phase transition begins, 
S{k,t) develops a peak at low k. As the phase transition proceeds, this peak grows in 
amplitude and redshifts until at late times it completely dominates the infrared portion of 
the spectrum. 


III. DEFECT DENSITY 

We have argued that the dynamics of a phase transition of a quantum field may be 
approximated by the dynamics of the two-point function at least for times less than the 
spinodal time. Using the two-loop 2PI equations of motion we obtained a numerical solution 
for the evolution of the two-point function. To observe the formation of correlated domains 
and topological defects it is necessary to identify these structures from the form of the 
two-point function. 

The identihcation of topological defects from the underlying quantum dynamics of the 
two-point function is a complicated problem. Intuitively, the existence of well-defined cor¬ 
related domains of true vacuum at the completion of the phase transition is similar to the 
domains that form in a condensed matter system like a ferromagnet. However, our system is 
described by a quantum field theory and the existence of a classical configuration of domains 
requires a quantum-to-classical transition, of which decoherence is an essential condition. 
The extraction of a positive-definite probability distribution, P[‘F], from the density matrix 
requires a quantitative treatment of the quantum-to-classical transition. Classical defect 
solutions would then appear in held conhgurations drawn from P[‘F] ||12|| . 
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In this research our focus is not on the extraction of classical defects from the quantum 
system but instead on the long wavelength modes that determine the size of correlated 
domains and, therefore, the average defect separation. The long wavelength modes interact 
with an environment of short wavelength and thermal fluctuations that destroys quantum 
coherence among the long wavelength modes and results in a hnite correlation length for 
the system [Q. While quantum fluctuations are indeed important to the description of 
phenomena within the scale of one domain, beyond this scale we may hope to treat the modes 
as classical. With this abbreviated estimate of decoherence we show that the existence and 
size of correlated domains is indicated by an infrared peak in the power spectrum of the 
equal-time momentum-space two-point function. 

For a free held theory quenched into the unstable region by an instantaneous change in 
the sign of the square of the mass at t = to h is possible to obtain analytic expressions for 
the equal-time momentum space two-point function After the quench, the momentum 
space structure function k‘^G{k, t) has a strong maximum at the value 


kr, 


2t ■ 


(3,1) 


The real space Fourier transform of this function (normalized to unity at 


= 0) is 


G( 




m ( I r — r’ ) 


/K2 


~ e 


(3.2) 


This form of the equal-time correlation function is common in condensed matter systems 
H- From the equal-time real space correlation function we can identify the correlation 


length 


m 


8t 


irif kr, 


(3,3) 


In analogy with the free held theory, we assume the domain size is proportional to the size 
of the maximum of the momentum space structure function. 

Even when the density matrix of the quantum held has decohered into a positive-dehnite 
probability distribution for an ensemble of classical held conhgurations, the identihcation 
of topological defect structures is a complicated task. In a classical held conhguration and 
when defects are well-formed (so that the width of the defect is much smaller than the 
typical defect spacing), it is possible to identify zeros of the classical held conhgurations 
with topological defects. In a model with a global 0{n) symmetry in n spatial dimensions, 
and when the held probability distribution is Gaussian, a formula for the ensemble average 
density of held zeros was given by Halperin and derived explicitly by Liu and Mazenko 


m-- 


Ci = 



(3.4) 
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For a single scalar field in 3 spatial dimensions Eq. ( |3.4|) with n = 1 is valid as the ensemble- 
averaged density of zeros along a one-dimensional section of the field. The validity of the 
Gaussian approximation is further discussed in . 


A zero of the classical held conhguration does not uniquely identify a topological defect. 
Thermal huctuations give a large number of zeros of the held conhguration on small scales 
^ . The zeros of the held conhguration which are not associated with diherent defects lead 
to an ultraviolet momentum divergence in the expression for the zero density, Eq. (0. 
To make physical sense of the zero formula, a coarse graining of the held conhguration is 
needed. The spacing of defects and the size of correlated domains is determined by long 
wavelength excitations and if zeros of the held conhguration are to accurately count defects, 
it is necessary to apply a coarse-graining that removes the short-wavelength excitations 
caused by quantum and thermal huctuations. To ehect this coarse-graining in Eq. ( |3.4|) we 
impose a spatial momentum cut-oh at the upper edge of the unstable momentum band. 
Only the unstable modes will contribute to the domain size. 

In principle it is possible that the length scales set by kmaxif) and by pit) are diherent. 
When the structure factor kfG{k, t) is very strongly peaked about kmax we can approximate 
the unstable portion of the spectrum as 


PG{k,t) ^ 5{k- kmax) 


the zero density for a single scalar held in 3 dimensions is then 




vr y/ k‘^G{k, t) 


(3.5) 


(3,6) 


For the parameters used in the simulation and at late times, the scales set by the zero density 
and by the domain size diher only by a numerical factor of order unity. Using the linear 
model introduced in Sec. V, evidence that the growth of domains is characterized by one 
scale is shown in Fig. 5. Additional evidence for a “one-scale” model was given in in 
the context of the dynamics of a classical held phase transition in 1 -|- 1 dimensions. 

For very weak coupling not considered here and when a peak is no longer evident in 
S{k,t), the one-scale approximation will fail and it is possible in principle for the defect 
density to depart signihcantly from one per correlation volume . 


IV. CRITICAL SCALING 

The evolution of the two-point function and the extraction of correlated domains and 
topological defects allows a hrst-principles analysis of the mechanisms important for topo¬ 
logical defect formation. In light of the “freeze-out” scenario it is interesting to compare 
the initial size of correlated domains for scale factors ar{t) with diherent expansion rates r. 
To quantify the dependence of the initial size of correlated domains on the quench rate r of 
the phase transition, we compare the size of domains for diherent values of the parameter 
r. As discussed above, the average size of a domain is proportional to the maximum in the 
infrared portion of S{k,t). We compare the maximum kmax for diherent values of r and at 
two distinct sets of times during the phase transition. The hrst set of domains is measured 
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when the square of the effective mass reaches a minimum value. This provides an early-time 
measure of the size of domains. The results are shown in Fig. |^. The second set of domains 
is measured when the square of the effective mass reaches a local maximum value after the 
phase transition. This provides a late-time measure of the size of domains. The results are 
shown in Fig. ^ and Fig |^. For slow quench rates (r > 1.0), at both early and late times, 
the dependence of kmax on the quench rate r is well approximated by a power law, 

kmax{r) ~ (4.1) 

For fast quench rates (r < 1.0), the dependence of kmax on the quench rate r is well 
approximated by a power law. 


kmax{r) ~ r (4.2) 

The origin of these exponents and the difference between slow and fast quenches is discussed 
in the following section. 


A. The Preeze-Out Scenario 


In Zurek’s scenario the “frozen” correlation length and therefore the initial size of corre¬ 
lated domains scale as a power law of the quench time r as 




7" i+fi 


(4.3) 


where p and u are the equilibrium critical exponents for the correlation length and relaxation 
time respectively. The correlation length and relaxation time are identihed as, respectively, 
the length and time scales that characterize the equilibrium behavior of the propagator near 
the critical temperature |^. Under the scaling hypothesis, the equal-time propagator is 
written 


Gk = k-^^-^^F[h{e)k]. 


(4.4) 


Here e is the reduced temperature, Eq. (pT^) , p is a critical exponent (not to be confused 
with conformal time) and F is a dimensionless function. The hnite-temperature equilibrium 
correlation length is 


^(e) = h{e) ~ e (4.5) 

Similarly the two-time propagator is 

G^^k=u;‘^F[h\e)u;,h{e)k] (4.6) 

so that the relaxation time 

T h' e~^. (4.7) 

In the Hartree-Fock approximation the equilibrium propagator is 
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m 


,/F + miff 

2 coth(-^^^--/?) 


(4.8) 


eff 


where the effective mass is given by the equilibrium value Eq. (|2.21|) . Near the critical point 




( 4 . 9 ) 


where A is a constant independent of temperature. The scaling behavior of the propagator 
is as ^ and and therefore the critical exponents for the theory are 




( 4 . 10 ) 


The critical exponents derived for this quantum field theory in 4 spacetime dimensions are 
the same as the critical exponents of those of a classical Ising model in 4 space dimensions 
since, after the Euclidean continuation t —> if3, they are in the same universality class. 
The finite size (3 of the Euclideanized time dimension can alter the critical behavior of the 
quantum held . At high temperature or near a second-order phase transition equilibrium 
quantum systems in d spacetime dimensions undergo dimensional reduction and behave with 
the critical exponents of classical systems in d — 1 space dimensions |4^. The dimensional 
crossover can be parameterized by 


z = ^tT, 


( 4 . 11 ) 


where is the equilibrium correlation length and T is the temperature. The effective critical 
exponents are functions of z so that z = 0 corresponds to the quantum case and z = oo 
corresponds to the classical case. Intuitively, dimensional crossover occurs in a second- 
order phase transition when the correlation length is much larger than the Euclidean 
compactified dimension (3 = 1/T. The quantitative effect of dimensional reduction on the 
critical exponents has been calculated in detail in where plots of Ueffi^) and p,eff{z) may 
be found. In the extreme dimensionally-reduced case, when z ^ oo, the critical exponents 
change to 


H = 1/ ze 0.64. (4-12) 

However, in our model, the dynamics of the system and in particular the “freeze-out” of the 
correlation length prevent z from becoming too large and the actual critical exponents he 
somewhere between the four dimensional case Eq. ( |4.1CI|) and the three dimensional case Eq. 
( |4.12| ). The “freeze-out” prediction for the initial scale of the density of topological defects, 
Eq. (^4.3|) , is relatively insensitive to this change in critical exponents and varies from 

(4.13) 


in the four dimensional case to 


e/ 


reeze 


m.38 
~ r 


( 4 . 14 ) 
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in the three dimensional case. Our observation of the power law exponent of 0.35 in the 
numerical simulations for slow quenches, r > 1.0, (Fig. 3 and Fig. 4) is consistent with 
these cases. 

The calculation of the equilibrium critical exponent for the relaxation time depends on 
the dynamics of the mode functions for the the low k modes near the critical point and in 
particular whether the dynamics is overdamped or underdamped. Near the critical point, 
mlj-f ~ 0, and the mode function equation is approximately 


(P a{tc) d \ 

a(tc) dt a^(tc)) 


fk{t) = 0. 


(4.15) 


The critical time tc can be estimated as the time when the equilibrium effective mass, Eq. 
( ^(^ , goes to zero. 


2 


(4.16) 


and the value of the scale factor at tr is 


a{tc) = 


tc + T 


T 


1.3. 


(4.17) 


The mode function equation near the critical point Eq. ([4.15|) is that of a damped harmonic 
oscillator with natural frequency (Uq = and damping constant F where 


The modes for which (Uq > F or 


F = 3H{tc) ^ -. 


k > — 

T 


(4.18) 


(4.19) 


are underdamped. The range of wavenumbers k for the modes that determinex the size of 
correlated domains also depends on r. Eq. ([4.19|) is a condition on r such that when 


r > r* 


(4.20) 


the modes responsible for domain formation are underdamped and r* is to be determined. 
For the underdamped case we expect kmax ~ r“3. Using kmax as a representative wavevector 
of the modes responsible for domains, Eq. ( [1.19|) implies the underdamped condition 


r* 


1 . 


(4.21) 


Slow quenches with r > 1.0 are therefore underdamped. Fast quenches are overdamped and 
the dynamics of the mode functions will be dominated by the hrst time derivative. In the 
overdamped case, the critical exponent, /r, now assumes the non-relativistic value 


l^overdamped ‘^k'underdamped- 


(4.22) 
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Due to the uncertainties of dimensional reduction, in the overdamped case, the initial scale 
of the density of topological defects ranges from 

^free.e ~ (4.23) 

for the four dimensional critical exponents to 

^freeze ~ (4.24) 

in the three dimensional case. This is in good agreement with the power law exponent of 
0.28 measured for overdamped, fast quenches (r < 1.0) shown in Fig. 

In both the overdamped and underdamped cases, the prediction of the “freeze-out” 
proposal for the power law exponent of the scaling of the initial defect density with the 
quench rate appears to be in excellent agreement with the numerical simulations. 


B. Early-Time Linear Approximation 


If the initial size of correlated domains and the initial density of topological defects are 
determined by processes occurring very near the critical point, as claimed in the “freeze-out” 
proposal, then the power law scaling of domains with the quench parameter r observed in 
the simulation of the Hartree-Fock mode function Eq. (|2.19| ) will also appear in a linear 
approximation. A linear amplitude approximation to the full mode function equation is valid 
for slow quenches and only near the onset of the spinodal instability, before the unstable 
modes have grown appreciably and back reaction is important. Consider the linear equation 



\d^ dd 2 I 

dt^ a dt 

Fk = 0 (4.25) 

where 

K^ = e + kl 

(4.26) 

and fco is related to the initial temperature 



^0 - 24^0 

(4.27) 

In conformal time, we have 



rf = At {t + t) 

(4.28) 

with conformal modes, fk, dehned by 



> 

III 

(4.29) 


Eq. (|4.25|) may be rewritten in conformal time as 

drj'^ 


- 


2 2 
m Tj 

4r^ 


/a: = 0. 


(4.30) 
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This equation can be solved exactly in terms of parabolic cylinder functions. However, a 
simpler solution is obtained for times near the critical point. We further approximate it by 




mK 

T 


iv - Vk) 


fk 


0 


Vk 


2Kt 

m 


(4.31) 

(4.32) 


The properly normalized solution with vacuum boundary conditions before the instability 
is 


fk = i 



2 

mKx^ 

3 \ 

T 


(4.33) 


where 


x = v -Vk 


(4.34) 


and i^i /3 is a modihed Bessel function. 

We now use the analytic solution Eq. (|4.33|) to examine the dependence of the position 
of the peak in the Fourier space structure factor k‘^G{k,t) on the quench parameter r. The 
location of this peak kmaxif) redshifts throughout the phase transition, moving towards lower 
momentum as the domains coarsen. In order to compare domains formed in the linear model 
with those of the Hartree evolution we compare the position of the peaks at very early times 
when a peak is hrst identihable in S{k,t). Plots of S{k,t) with the analytic modes were 
made with Mathematica. The results are shown in Fig. |^. The scaling of domains with 
the quench parameter r with the same power law exponent as seen in the full numerical 
simulation for slow quenches is evident. 


V. DISCUSSION AND CONCLUSION 


Using the Hartree-Fock approximation to the equations of motion for the two-point 
function of a quantum scalar held undergoing a phase transition in a 3-1-1 dimensional 
spatially hat, radiation-dominated FRW Universe we have shown that the size of correlated 
domains, measured as the maximum of the peak of the infrared part of the spectrum of 
k'^G{k,t), scales as a power of the quench rate r as 




omains 


.0.35 


(5.1) 


for slow, underdamped quenches (r > 1.0) and 


id 


omains 


.0.28 


r 


(5.2) 


for fast, overdamped quenches (r < 1.0). The observed power law scaling of correlated 
domains is quantitatively consistent with the “freeze-out” hypothesis. In both overdamped 
and underdamped cases, the value of the power law exponent extracted from the Hartree- 
Fock evolution of the two-point function is in good agreement with the value calculated in the 
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“freeze-out” scenario using the critical scaling exponents for a theory in the Hartree-Fock 
approximation. 

To further explore the behavior of the quantum system near the critical point we in¬ 
troduced an approximate linear model valid for slow quenches and short times after the 
onset of spinodal instability. In this linear model, the size of correlated domains scales with 
the quench rate r to the same power as observed in the underdamped numerical simula¬ 
tions and predicted by the “freeze-out” proposal. This provides analytical evidence for the 
“freeze-out” hypothesis. A similar analytic model was introduced in the context of a 1-1-1 
dimensional classical condensed matter system in [Q 

In contrast to previous approaches to the problem of defect formation in quantum held 
theory, we have allowed for the process of back reaction, which permits the quantum held 
to exit the region of spinodal instability. The power law exponent is the same whether it is 
measured at very early times with the analytic linear model (Fig. |^) or, in the numerical 
simulations, at early times in the middle of the spinodal region at the minimum value of 
Mgjj (Fig. H) or, at late times in the stable region at the maximum value of (Fig. ||). 
This is consistent with the idea that the relevant processes for the growth of domains occur 
near the critical point. However, back reaction is essential to the freezing in of the value 
of the power law exponent and to the “freeze-out” hypothesis. Although the power law is 
accurately recorded in the early time analytical model, if back reaction is ignored and the 
linear model is (incorrectly) extrapolated to late times, the predicted scaling exponent is 
different 


kmax^T 2. (5.3) 

The results derived from the numerical simulation of the mode function equation with 
back reaction and the analytical solution around the critical point are evidence supporting 
the hrst verihcation of the “freeze-out” scenario in 3 spatial dimensions and in a realistic 
system relevant to the early Universe. It is at hrst surprising that arguments based on 
classical equilibrium critical scaling apply in the context of dynamical quantum held theory. 
A possible explanation rests with the high temperature initial conditions and the properties 
of classical theory in FRW spacetime. Since the initial temperature is much higher than 
the initial ehective mass 


To >> me//(0) (5.4) 

the held is approximately conformally invariant. In an FRW Universe, a conformally invari¬ 
ant held initially in thermal equilibrium will remain in equilibrium at a redshifted temper¬ 
ature 


m = ( 5 . 5 ) 

a[t) 

A well-known example is the redshifting of the blackbody spectrum of the cosmic microwave 
background radiation as the Universe expands. To illustrate the approximate conformal 
invariance in our model we compare the dynamical ehective mass from the simulations with 
the ehective mass of a theory with the equilibrium redshifted temperature of Eq. (|5.5|) 
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Meq = -m 
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-'o 
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The results are shown in Fig. Only after the onset of the spinodal instability do the 
equilibrium and dynamical effective mass differ signihcantly. The approximate conformal 
invariance of the theory means that the finite-temperature effective potential and the critical 
behavior derived from it are approximately valid until times near the onset of the spinodal 
instability. The high temperature initial conditions also offer a possible explanation for the 
observed classical behavior of the system. At such high temperature, thermal fluctuations 
dominate over quantum vacuum fluctuations. 

The agreement in the value of the power law exponent between the microscopic evolu¬ 
tion equations of the quantum field theory and phenomenological critical scaling supports 
the contention that quantum critical systems in the early Universe share, in certain circum¬ 
stances, many of the properties of their classical counterparts. It would be very interesting 
to further explore these connections. To do so, however, it is necessary to go beyond the 
strong coupling and high temperature conditions used in this research. 

The results derived in this paper represent a preliminary step in a hrst-principles ap¬ 
proach to the calculation of the topological defect density immediately following the comple¬ 
tion of a second-order phase transition in the early Universe. Our work uses a microscopic 
quantum field theory and incorporates realistic initial conditions and a physical quench 
mechanism. However, the domain wall defects formed in this model are inconsistent with 
cosmological observations. It is not difficult in principle to apply the methods used in this 
paper to more cosmologically realistic theories such as cosmic string models. The analytical 
and numerical evidence for the power law scaling of the domains is expected to hold in a 
more realistic model. 

The formalism of the 2PI-CTP effective action and related techniques can be used to 
probe more general questions related to the physical aspects of quantum critical dynamics 
m- To observe and isolate quantum processes it is necessary to relax the assumption 


of high temperature initial conditions, allowing quantum vacuum fluctuations to dominate 
over thermal fluctuations. Also important are the detailed properties of the system-bath 
interaction. A system-bath interaction such as proposed in could address such issues as 
critical slowing down and the role of dissipation and noise. Such work is in progress. 
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FIG. 1. Plot of the square of the (renormalized) effective mass M = vs. cosmological 

time t for quench parameter r = 1.0. 



k 


FIG. 2. Plot of the Fourier-space structure factor, S{k, t) = k'^G{k, t), vs. comoving momentum 
k at various times during the evolution for quench parameter r = 1.0. The bottom curve is a 
snapshot of S{k,t) at time t = 3.9. The middle curve is a snapshot at t = 5.4. The top cnrve is a 
snapshot at t = 6.9. 
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FIG. 3. Plot of logioC^max) VS. logio('7‘) in the case of slow, underdamped quenches (r > 1.0) 
for domains formed early in the phase transition, at the time when the square of the effective mass 
reaches a minimum value. Filled squares are measurements from the numerical simulations. The 
solid line is a plot of the best-fit linear function to the data, \ogiQ{kmax) = ~0-35 logio('^) “ 0.14. 



Log(tau) 

FIG. 4. Plot of \ogiQ{kmax) VS. logio('^) ™ of slow, underdamped quenches (r > 1.0) 

for domains formed late in the phase transition, at the time when the square of the effective mass 
reaches a local maximum value. Filled squares are measurements from the numerical simulations. 
The solid line is a plot of the best-fit linear function to the data, \ogiQ{kmax) = —0.35 log]^Q(r)—0.26. 
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FIG. 5. Plot of \ogiQ{km,ax) VS. log^o(''’) ™ the case of fast, overdamped quenches (r < 1.0) 
for domains formed late in the phase transition, at the time when the square of the effective mass 
reaches a local maximum value. Filled squares are measurements from the numerical simulations. 
The solid line is a plot of the best-fit linear function to the data, \ogiQ{kmax) = —0.28 log]^Q(r)—0.22. 
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FIG. 6. Plot of the ratio, R{t) = 


p{t) 


.it) 


vs. cosmological time t for quench parameter 


r = 10.0. The defect density p{t) was calculated using the analytic modes, Eq.( 4.33| ), and the 
Halperin-Mazenko-Liu formula, Eq. (^), with a cutoff at the maximum momentum of the unsta¬ 
ble band. The maximum kmax{t) in the structure function, S{k, t) = k‘^G{k, t), was also determined 
using the analytic modes, Eq.( 4.33| ). At late times this ratio approaches a constant that is inde¬ 
pendent of the quench parameter. 



Log{tau) 


FIG. 7. Plot of )ogiQ{kmax) VS. log]^Q(T) in the analytic model, Eq. ( 4.33 ). The maximum, 
kmaxt was determined as soon as a peak was evident in the structure function, S{k,t). Filled 
squares are measurements from plots of the analytic modes. The solid line is a plot of the best fit 
linear function to the data, \ogiQ{kmax) = ~0.341ogio(T) — 0.09. 


30 













FIG. 8. Plot of the square of the effective mass M = vs. cosmological time t in the 

equilibrium case (dashed line) where ~ —and the in nonequilibrium case (solid 
line) where w-g// is determined by the full numerical simulations. The slight difference in the 
equilibrium and nonequilibrium curves at the beginning of the evolution is due to the difference 
between the general Hartree-Fock effective mass given by Eq. ( 2.20| ) and the high temperature, 
small A limit given by Eq. (p.21). 
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